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ABSTRACT 



Aims. We investigate a new approach to the detection of companions to extrasolar planets beyond the transit method. We discuss the 
possibility of the existence of "binary planets". 

Methods. We develop a method based on the imaging of a planet-companion as an unresolved system (but resolved from its parent 
star). It makes use of planet-companion "mutual phenomena", namely mutual transits and mutual shadows. 
Results. We show that companions can be detected and their radius measured down to lunar sizes. 

Key words, (stars:) planetary systems 



1. Introduction 

The search for satellites of extrasolar planets is relevant to the 
understanding of the evolution of planetary systems and to the 
perspective of their habitability. Their occurrence in other plane- 
tary systems is very likely since in the Solar System, 7 of the 
9 "planets" have from 1 to several tens of satellites and nu- 
merical simulations of planet formation show that such satel- 
lites should be common (Ida et al. 11997) . Several methods have 
been proposed for their detection: Sartoretti & Schneider ( 1999 ) 
proposed detection by transits; Brown et al. (2001) applied this 
approach to constrain companions to HD 209458 b; Han & 
Han (2002 ) and Bennett & Rhie (2002 ) discussed the possibility 
of detecting companions via microlensing; Williams & Knacke 
(120041 ) have shown that Earth-like satellites of giant planets in 
the habitable zone would be detectable by spectroscopy in the 
CH4 hole of the planet spectrum in the 1.5- 4fim band. Here we 
investigate new ways to detect planet satellites based on planet 
wobble and photometry in direct imaging. First, we discuss the 
possibility of binary planets and their impact on planet charac- 
terization. 

In addition to satellites defined as companions with masses 
and sizes significantly smaller than their parent planets, we sug- 
gest that there may exist binary planets with comparable masses 
and sizes. Such binary planets do not exist in the Solar System. 
However, the detection of the first 180 other planetary systems 
has shown much diversity and differences with the Solar System 
in their characteristics and provided surprises such as very small 
and highly eccentric orbits. These unexpected findings open the 
possibility of new configurations such as binary planets. In ad- 
dition, binarity is not an exception in objects conceptually as- 
sociated with planets, such as asteroids (Pravec et al. 2006), 
trans-Neptunian objects (Kern & Elliot 2006) or brown dwarfs 
(Stassun, Matthieu & Valenti 120061 ). In the same way, refined 
simulations show the formation of vortices close to planets (Kley 
2003 ) which could in turn form a massive companion to the 
planet. 

An interesting counterargument has recently been put for- 
ward (Canup & Ward 2006 ) predicting that the mass ratio of the 
companion to its host planet cannot be larger than about 7 10~ 4 . 



It will be interesting to see in the coming decade if this predic- 
tion is confirmed by observations. We note that this prediction 
does not hold for non giant Solar System bodies like the Earth 
and Pluto. 

The presence of companions to planets would have a signif- 
icant impact on their characterization. 

- Planet mass 

For binary planets the radial velocity method of detection 
(RV) gives only the sum M — M\ + M2 of their (minimum) 
masses, leading to a false assignment of the mass of indi- 
vidual objects. In particular if a companion detected by RV 
has a (minimum) mass larger than the standard planet mass 
upper limit (~ 13 Jup. mass), it would be inappropriately 
discarded as being a brown dwarf and not a planet. The same 
consideration holds for the astrometric detection of planets. 

- Radius and albedo 

The future detection of planets by direct imaging will give, 
in the case of a planet-companion system, only the sum 
F = F\ + F2 of their fluxes: 

(^Refl)i,2 = A\,2 x ^12 (modulated by an orbital phase) for 
the reflected flux and 

(^Th)i,2 = R\ 2 x T* 2 (constant along the orbit) for the ther- 
mal emission. 

The albedo A of a cold planet cannot be larger than 1; 
its radius cannot be larger than ~ 1.1 Rj up (Guillot 2005 ). 
Therefore the normalized flux of a planet cannot be larger 
than the maximum normalized flux: 



Fmax = 1.2 R Jnp and F max 



= r 4 i.2 r 2 

Jup 



for the reflected and thermal flux respectively. 
If the observed normalized flux of the planet is larger than 
the maximum possible flux F max , either there is something 
wrong with its structure or it is binary. An example of an odd 
structure is given by a planet surrounded by rings, since then 
there is an additional contribution to the planet flux coming 
from the planetary ring (Schneider 2003). Their presence can 
be detected in reflected light by a non-Keplerian modulation 
of their reflected flux along the orbit (Arnold & Schneider 
120041 Barnes & Fortney 2004); in the thermal regime there 
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is no such orbital modulation and there is a priori no way to 
disentangle a binary planet from a ringed planet. 
If the observed planet normalized flux is smaller than F max 
and if it is in reality a planet-companion system, the assump- 
tion that it is a single planet would lead to incorrect assign- 
ments to the albedo and/or radius of its components. 
- Planet spectrum 

If there is a satellite with a significant contribution to the 
planet flux, the global planet + satellite spectrum may lead 
to misinterpretations if the spectrum is attributed to a single 
object (leading to incorrect atmosphere models). An impor- 
tant example is provided by an icy satellite of a telluric planet 
close to the habitable zone. The planet may undergo a green- 
house effect leading to an atmosphere with water and C02, 
while the satellite can be icy. In that case, photodissociation 
of the satellite's ice can lead to ozone synthesis (Teolis et al. 
2006); the global spectrum would mimic an Earth spectrum. 

The case of planets detected by transits is different. Their 
possible binarity would be easily inferred from the peculiar 
shape of the transit lightcurve or from their timing revealing 
non periodic transits (Sartoretti & Schneider 1999). For the 
9 planets detected by transit by 22 Aug. 2005 (HD 209458 
b, HD 149026 b, TRES-1, OGLE-TR-10, OGLE-TR-111, 
OGLE-TR-1 13, OGLE-TR-132, OGLE-TR-56, HD 189733; see 
exoplanet . eu), none of them presents any sign of a moon or of 
binarity. But their orbital radius is so small (a < 0.05 AU) that the 
Hill radius a(M p i/3M„) 1/3 (Valtonen l20061 l inside which a stable 
orbit can survive around the planet would make these planets 
merge together (for a study of the stability of satellites around 
giant extrasolar planets, see Barnes & O'Brien 2002). Binarity 
among "hot Jupiters" is therefore very unlikely. Only planets at 
orbital distances larger than ~ 0.1 AU can, from a dynamical sta- 
bility point of view, be binary. But the probability that they make 
a transit drops as I /a and only large surveys like the CoRoT and 
Kepler space missions will detect them. Finally, the detection 
by transits of a binary planet lacks generality since it requires a 
transit to occur. 

Thus, the detection of planets by RV, astrometry and direct 
imaging miss their possible binarity and could lead to misas- 
signments of their mass, albedo and radius. The most radical 
method to detected the binarity of a planet would consist of re- 
solving the planet-companion system in high angular resolution 
direct imaging. But, for a typical Jupiter/Ganymede system at 
5pc, the baseline B required to resolve the system at 2/rni would 
be B = 1.2 x 2yum x 5pc/10 6 km = 360m. In addition, a very 
high contrast of at least 10 6 is required. Such a baseline and high 
contrast only will be achieved in the future, here we investigate 
less difficult methods. 

2. New approaches to planet-companion system 
detection 

We consider the detection of a planet-companion system by di- 
rect imaging in which the system is unresolved and thus appears 
as a single point, supposed to be detached from its parent star. 
Gaidos et al. (2006) have considered the thermal emission ap- 
proach. Here we concentrate on the reflected light approach of 
planet imaging. 

2.1. Planet wobble 

Let F p i = 7rA p i/?pj and F c = nA c Rl be the flux of the two compo- 
nents of the planet-companion system. The photocenter will then 



make a wobble around the center of mass with an amplitude Act 
given by 



Aa = 



M p i + M c 



F pl M c - F C M, 



f pi + F c 



(1) 



For a Saturn-Titan system, the mass and flux ratios are 
MJM pi = 0.01 and FJF pi = 0.001; then at 5 pc, the ampli- 
tude of the angular wobble is 15 microarcsec. For a Jupiter- 
Saturn binary system with the same separation M c /M p i = 0.3 
and F c /F p \ = 1; the angular wobble amplitude is then 1 mas. 
For a perfect twin planet (same mass, radius and albedo), the 
photocenter remains fixed along its revolution and the binarity is 
undetectable by this approach. 

These numbers should be compared to the accuracy 66 on 
the planet position from ground and space observations. This 
accuracy is given, at best, at a wavelength A, for a bandwidth W 
and a telescope of diameter d by: 



66 = 



A 

'N P h ~ d 



(2) 



where N ph = W W^ Am (d/1 cm) z (W/1A) (T exp /1 sec) is 
the number of detected photons emitted by the planet. For a 
m = 25 planet, 66 = 15 //as (resp. 1 //as) in a 10 h exposure with 
W = 100 nm and d = 8 m (resp. d = 30 m) at A = 1 /an. 

Of course, any intrinsic flux variation of the planet and/or 
companion would complicate the detection of the signal. 

The planet-companion system can be seen as a spectroscopic 
binary. Each component has a radial velocity variation with an 
amplitude V u = [M 2 ,\l{M l + M 2 )] (G (Mi + M 2 )/a c ) 1/2 . 

For a Saturn-Titan system, the exo-Saturn velocity is 
2 m/sec. For a Jupiter-Saturn binary system with the same sepa- 
ration the exo-Jupiter velocity is 10 km/sec. 

The feasibility of this wobble approach has been addressed 
in the context of Extremely Large Telescopes by Ardeberg 
(2005). 

2.2. Mutual photometric phenomena 

Let us consider a restricted three body problem consisting of 
a star, a planet orbiting the star and a companion orbiting the 
planet. In our problem, the companion (with radius R c and mass 
M c ) and the planet (with radius R p \ and mass M p i) orbit each 
other in a circular orbit of semimajor axis a c ; the center of mass 
of this system orbits the star (with radius R* and mass M„) with 
a semimajor axis of a p \ (see Fig. [J}. 

The only source of light is the star, which shines with con- 
stant luminosity P (following the notation of Lester et al. |1979) . 
The flux (density) F, arriving at the planet and the companion is 
constant and uniform. The flux reflected by these bodies is not 
resolved, but it is detached from the light coming from the star. 
This reflected flux depends on the physical characteristics of the 
reflecting surfaces and the geometrical configuration of the sys- 
tem. 

The theory that describes the radiation by a planet illumi- 
nated by its star is well described in Lester et al. (1 19791 1 and 
Fairbairn (2002) for example. Here we will suppose that both 
the planet and the satellite are Lambertian spheres. Scattering 
atmospheres may differ from this model but ours is still a good 
approximation. 

In these circumstances, the flux reflected by the planet (see 
Eq. |4]i is proportional to: 
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alpha = 
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alpha = pi/2 
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alpha = pi/3 
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alpha = 2pi/3 



Fig. 2. Illumination of the planet as seen from the observer for 
different phase angles. 



PHASE 1: 
companions's shadow on the planet 



Fig. 1. Schematic view of the system (not to scale). The compan- 
ion's orbit around the planet and planet's orbit around the star do 
not need to be coplanar. The white surfaces of the planet and the 
companion are illuminated, while black surfaces are in darkness. 



PHASE 2: 

companion transiting in front of the planet 
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- its geometrical albedo p; 

- its surface (that is, the square of the radius R); 

- the inverse of the square of its distance to the star a; 

- the phase law (a). 

The phase law is an integral which takes into account the ge- 
ometry of the system. It depends on the phase angle a, which is 
the angle between the incident direction of the radiation and the 
direction to the observer, as seen from the center of the sphere 
(see Fig. [2j». The integral has the following form: 



d6d<p sin (&) cos(0) cos(0 + a). 



(3) 



With these assumptions, the flux reflected by the planet as 
seen from Earth is: 



R 2 

F p i=p-f <F(a)F» 
fl pi 



(4) 



The limits of integration of Eq. (01 are 8 e [0, n] and (f> e 
[-n/2, jt/2 + a]; so it takes the form: 



sin(a) + (n - a) cos(a') 



(5) 



For further details on how to arrive at these expressions, 
please see Lester et al. (119791 1. 

In general, the total flux that we measure is the addition of 
the individual fluxes coming from each body. However, if the 
orbit of the satellite lies on the line of sight of the observer, the 
orbital revolution of the satellite around the planet has 4 remark- 
able phases (see Fig. [3): 

- Phase 1 : companion's shadow on the planet. 

- Phase 2: companion transiting in front of the planet. 

- Phase 3: companion eclipsed by the planet. 



PHASE 3: 
companion eclipsed by the planet 



PHASE 4: 
companion occulted by the planet 





Fig. 3. Mutual phenomena during the revolution of a compan- 
ion around its parent planet; seen from the point of view of the 
observer; illumination comes from the left. 

- Phase 4: companion passing behind the planet and being oc- 
culted by it. 

In any of these cases we see a flux decrease because part of 
the reflecting surface is in darkness. 

To calculate this flux decrease we consider the surface in 
darkness and we integrate Eq. © within the limits of this sur- 
face. We have calculated this flux numerically and we present it 
in Figs. |4] to |7] for different configurations of the system. In the 
next subsections we will detail these calculations. The signal in 
the lightcurves is not periodic any more (as it was in the transits 
of a single planet in front of a star); the position, shape and depth 
of the flux decrease strongly depend on the relative positions of 
all the bodies involved (star, planet, companion and observer). 

Mutual transits 

One of the bodies is transiting in front of the other, occulting 
its surface. This surface is the intersection of a cylinder with a 
sphere: a classic problem of surface geometry. However, the ex- 
pression resulting from integration © within the limits imposed 
by the intersection is not analytic. In the next two points we will 
approximate that integral for the two possible transits, in order 
of difficulty. 
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companion eclipsed 



companion occulted 
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companion eclipsed 



companion occulted 
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total flux (planet + satellite) 



total flux (planet + satellite) 




Fig. 4. Lightcurve for a planet with a companion; the phase angle 
is a = 90°; R c /R p \ = 1/3; the trend in the out-of-event part of 
the curve is due to the variation in phase of the planet, as it orbits 
the star. 



companion eclipsed companion occulted 



total flux (planet + satellite) 




Fig. 5. Lightcurve for a planet with a companion; the phase angle 
is a — 90°; R c /R p i = 1/3; the orbit of the companion is inclined 
85° with respect to the line of sight. 

planet occulted companion eclipsed companion occulted planet eclipsed 
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total flux (planet + satellite) 
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Fig. 6. Lightcurve for a planet with a companion; the phase angle 
isa= l35°;RJR pl = 1/3 
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Fig. 7. Lightcurves for a planet with a companion; the phase an- 
gle is a = 45°; fl c //? p i = 1/3 



The geometric probability of transits is R p \/a c « 10%. 

- Phase 4: companion occulted by the planet 

The companion is supposed to be smaller than the planet, so 
except in the ingress and egress phases, the companion is com- 
pletely occultecQ. Using Eq. ©, the corresponding flux decrease 
is: 



AFt 



PcRl 



F T F pi + F c p pl R 2 pl +p c Rl 

The duration of this transit is: 
PcR P i 



AT, 



- Phase 2: companion transiting in front of the planet 

We can consider the maximal occultation: this is, the maxi- 
mum (or minimum, depending on the reader's choice) of the flux 
decrease. To do so, we suppose that the surface occultecQ has a 
constant brightness (which is not true) and we approximate this 
area by the surface of a disk with the companion's radius. The 
corresponding flux decrease is: 

AF T PpiJ?^|cosar/2 
~ * (PpiR^+PcRlma) 

and the event duration is: 

P c R p \(l + cos a) 
Ar 2 = . 

2na c 

Mutual shadows 

1 if the plane of the orbit of the satellite is inclined, this may not be 
true (see Fig. [5}; in that case we will have to calculate which is exactly 
the occulted surface of the companion. 

2 regarding Eq. {3]( we can easily verify that the coordinates of the 
brightest point of the surface are 6 = n/2 and <f> = -a/2. 
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In phases 1 and 3 one body eclipses the other. In eclipses 
we have to differentiate the umbra (defined as the region of an 
eclipse where the light coming from the primary source is com- 
pletely excluded) from the penumbra (the region where only part 
of the light from the primary source is excluded). However, in the 
case of planets orbiting at distances similar to those of Jupiter or 
Saturn, the fraction of time in which the companion is only in 
penumbra compared with the fraction of time in which it is in 
umbra is of the order of 10 4 (whereas for the Earth-Moon sys- 
tem it is 0.7); so we decided to neglect ifl- 

In this case the surface of interest is the intersection of a 
cone with a sphere. In the next two points we will approximate 
the integral @ for the eclipses. 

If the companion's orbital plane is close to the planet ecliptic, 
the geometric probability of mutual shadows is close to 100%. 

- Phase 3: companion eclipsed by the planet 

The companion disappears completely in the shadow of the 
parent planet. The relative flux variation is then: 



AF T 



PcRl 



F T F p i + F c p pI #2 + pc R 2 
with a duration of: 



A7\ = 



is the radius of the perpendicular section of the planet's 
shadow cone at the position of the companion: 



R* — R 



a c =R p Al- 



R*-R 



pi a c 



flpi ' \ R p \ flpi 

Phase 1: satellite's shadow on the planet 



In this case the companion projects a shadow on the surface 
of the planet. To calculate the flux decrease we proceed as in 
phase 2 and we obtain: 

2 3 , 



AF T P P \ Rc J cos a l 2 



Ft ( Pp iR 2 pl +p c R 2 c )y( a ) 



where R' c is calculated in the same way as R' pl in phase 3. The 
duration of the passage is: 



Ar, 



P c Rp\(l +COSQ') 

2na c 



when the planet or companion surface is solid (such as during 
the lunar eclipses). The flux decrease is then identical to the re- 
flected light case. When the planet or its companion have a thick 
atmosphere, its thermal inertia and atmospheric circulation in- 
hibit the temperature drop and no flux decrease is seen. For thin 
atmospheres the situation is intermediate. The amount of thermal 
infrared flux decrease during shadows and eclipses thus provides 
a way to estimate the thermal inertia of the surface of the planet 
and its companion (Spencer |l9871 l. 



3. Feasibility and search strategy 

3. 1 . Photometric accuracy and detection limits 

For a companion to be detectable, the photometric accuracy must 
be better than the transit or shadow decrease in an exposure time 
shorter than, say, half the duration of the phenomenon. The du- 
rations are typically of the order of F comp /?p//7ra comp , i.e. 8 h 
for a Saturn-Titan system and 4 h for an Earth-Moon system. 
The photometric accuracy on the detection of a planet is not 
controlled by its own photon noise, but by the photon noise of 
the background consisting essentially, in real situations, of the 
speckles of the parent star halo. 

In current coronagraphic detection projects of planets in the 
visible, the speckle background is typically 100 times the planet 
signal, at an angular separation of 1 .5 to 2 A/D (D = telescope di- 
ameter) for a m = 5 star and coronagraphic rejection factor of the 
star of 50 000 on the stellar peak. This figure holds for Jupiter- 
like planets at 1 AU (planet to star flux ration 10~ 8 ) detected by a 
1.5 m telescope and for Earth-like planets at 1 AU (planet to star 
flux ratio 10 -10 ) detected by a 7.5 m TPF-C. 

For a companion with one third the planet radius the relative 
flux decrease during shadows and transits is ~ 10%, i.e. 10~ 3 the 
speckle background. The detector must collect 10 7 speckle + 
planet photons to detect the companion with a SNR = 3 in the 
speckle background. For a Jupiter-sized planet the accumulated 
exposure time required is 30 h for am = 5 star with a 1 .5 m tele- 
scope. Since this exposure time exceeds the duration of the event 
(8 h in case of a Saturn/Titan-like system) by a factor of 4, the 
companion is detectable in a continuous planet monitoring only 
in cumulative exposures over 4 revolutions around the planet, 
i.e. after 2 months of continuous monitoring. This is for instance 
the type of monitoring planned in the 'Super-Earth Explorer' 
project (Schneider et al. 2006). For TPF-C, the collecting area 
is 20 times larger and correspondingly the exposure time to de- 
tect a 0.37?j up companion drops to 30h/20 = 1.5 h. Note that a 
continuous monitoring is also required for the detection of sur- 
face inhomogeneities of planets (Ford et al. 1200 11 1. 



The light curves presented here are free of noise, which ob- 
viously would make the detection of these events more difficult. 

Mutual phenomena in the thermal regime 

In the thermal regime, where one detects the thermal emis- 
sion of a planet and its companion, mutual transits occur in 
the same way as for reflected light and the corresponding light 
curves have the same shapes. For shadows and eclipses the situ- 
ation is different. The flux decrease is related to the decrease in 
temperature of the part of the planet and its companion not illu- 
minated by the parent star. This drop in temperature is sudden 



3 The effect of penumbra will be to soften the sharpness of the egress 
and ingress, which are very difficult to observe anyway. 



3.2. Geometric probability of shadows and transits and 
observation strategy 

The geometric probability of mutual transits is p - /? p i/a S at- For 
a Saturn/Titan-like configuration it is approximately 10%. The 
probability that a planet and its companion make mutual eclipses 
(with respect to the parent star) depends on the inclination of the 
orbital plane of the planet-companion system with respect to the 
orbit of the system around the parent star. If both objects have 
been formed by accretion in a protoplanetary disk, it it is likely 
that the two orbital planes are close to each other. Thus mutual 
shadows should occur at any position of the planet-companion 
orbit around the star. Consequently, the geometric probability of 
mutual shadows should be nearly 100%. 
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For mutual shadows, which have the maximum geometric 
probability of occurrence, the total fraction of time for which 
the event occurs is: 

Ar,+Ar 3 3 flpi 

— ~ ■ (6) 

" comp ft ^comp 

It is of the order of 10% for a Saturn-Titan system. It is thus 
necessary to have a duty cycle of at least 90% in order not to 
miss this event. The duration of this high duty cycle must be at 
least P C omp, i-e- 15 days for a Saturn-Titan system. In conclu- 
sion, the detection of companions by mutual shadows (the most 
probable event) requires continuous imaging of the planet for at 
least about two weeks. Such continuous imaging is possible only 
from space. 
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